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Abstract. The strong interaction shift e^ D and broadening P[ D in pionic deuterium have been remeasured with 
high statistics by means of the 7rD(3p - Is) X-ray transition using the cyclotron trap and a high-resolution crystal 
spectrometer. Preliminary results are e^f = (-2325 ± 31) meV (repulsive) and P^f (1 171 meV which yields 

for the 7rD scattering length a nD = [-(24.8 ± 0.3) + ( (6.3 )] • 10~ 3 m' 1 . From the imaginary part the threshold 

parameter for pion production is determined to be a = (252 * u ) pb. 



1 Introduction 

The deuterium experiment is embedded in a series of mea- 
surements |Q][2][3| aiming at a precision determination of 
the pion-nucleon scattering length comprising 

- pionic hydrogen (ttH), 

- pionic deuterium (ttD), 

- muonic hydrogen (juH), and 

- the characterisation of the Bragg crystal spectrometer 
using narrow X-ray lines from few-electron atoms pro- 
duced by means of an electron-cyclotron resonance trap 
(ECRIT). 

In this contribution, first results of the 7rD measurement are 
presented. 



2 Strong-interaction effects 

The complex pion-nucleus scattering length a„A of a pionic 
atom formed with a nucleus A(Z, N) is related in leading 
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order to the ls-state shift e\ s and width Fi, by the Deser 
formula |@) 

.T U 2aVc 4 

e\ s -i— = 7 a„ A + - (1) 

2 nc 

Ellipsis stand for corrections to the Deser formula because 
a K o is determined from a Coulomb bound state . 

In 7rH, the ground-state shift and broadening may be 
related to the nN isoscalar and isovector scattering lengths 
a + and a~ by 

ejf oc Orp-Hrp = a + +aT + ... 
rf s oc (a„- p ^o n ) 2 = 2( fl -) 2 + ... . (2) 

In this case, r oc a 2 because the hadronic broadening is 
due to a scattering process - the charge-exchange reaction 
n~p — > n n [6|. The corrections here have been calculated 
in various approaches and amount up to a few per cent of 
the leading contribution and include to some extent already 
electromagnetic and strong isospin-breaking terms [7,8,9 

isanna. 

In the case of 7rD, the real part 91 a„o of the scattering 
length, being proportional to the hadronic shift , maybe 
written by regarding the deuteron as a free proton and neu- 
tron as leading order and then applying additional correc- 
tions. One may write 

2a + + ... (3) 



EPJ Web of Conferences 



where ellipsis stand for the higher order terms from mul- 
tiple scattering (depending also on a + and aT), absorptive 
and possibly electromagnetic and strong isospin-breaking 
corrections. The second order correction is dominated by 
(cT) 2 being comparable in magnitude to the leading order 
term 11311141 because of the smallness of a + as required 
by chiral symmetry II15II161 . The nuclear structure is taken 
into account by folding with the deuteron wave function. 

The pionic deuterium shift provides a constraint on the 
jiN scattering lengths a + and cT . Its quality decisively de- 
pends on the precision of the data triple (e^ D ,e^ H ,r^ H ), for 
which a first set of precision data has been obtained in the 
last decade [ 17 18 19 20 |. Vice versa, starting from the el- 
ementary pion-nucleon processes and involving properly 
nuclear structure, multiple scattering, and absorptive phe- 
nomena, the hadronic s-level shift of the nD system must 
be calculable unambiguously ll2T1l22ll23ll24ll251l26ll27ll28l 

The forthcoming results for jtH from this measurement 
series will improve the constraint] 12 1, in particular by a 
better determination of the hadronic broadening in nH [3|. 
Finally, the accuracy for a + will be determined mainly by 
the poorly known low-energy constant f\ Ifl2l . 

Secondly, 7rD gives access to pion absorption and pro- 
duction AW <-> NNn at threshold 11321 . The hadronic broad- 
ening mainly is due to true absorption (n~d — > nn) be- 
ing the inverse and charge symmetric reaction of pion pro- 
duction {pp — > n + d) - in contrast to nH where the width is 
exclusively due to charge-exchange and radiative capture. 

The relative strength of true absorption to radiative cap- 
ture (n~d — > nny) was measured to S = — = 2.83 + 

v ' / nny 

0.04 1331 . and the branching ratio of internal pair conver- 
sion (n~d — > nne + e~) and charge exchange (ri~d — » nnn ) 
was found to be 0.7% J33 and 1.45 + 0.19 • 10~ 4 ED, re- 
spectively. Hence, the relative strength of the true absorp- 
tion channel to all other processes is obtained to S' = 
nn/(nny + nne + e~ + nnn°) = 2.76 ± 0.04, i. e., about 2/3 
of the hadronic width and with that of the imaginary part 
3 a n o is related to the process n~d —> nn, which can be 
linked with pion production. 

Pion production at low energies is usually parametrised 

by ma 

o-pp^ n+d = aClrj + /3C;y + ... (4) 

where 77 = p* n IM n is the reduced momentum of the pion in 
the nd rest frame. For 77 — > higher partial waves (J3, ...) 
vanish and only the threshold parameter a contributes rep- 
resenting pure s-wave production. The correction factors 
Ci take into account the Coulomb interaction. C\ is of the 
order of 30% [37,38 1 and an important source of uncer- 
tainty in the determination of a from cross-section data. 

Another approach to determine a is to exploit the nD 
ground state broadening where uncertainties stemming from 
Coulomb correction factors and normalisation of cross sec- 
tions are avoided. To derive the relation between a and 
3a„o as obtained from pionic deuterium, purely hadronic 
(non-experimental) cross sections <x are introduced to cir- 
cumvent the problem of diverging Coulomb cross section 
at threshold. Detailed balance relates pion production and 
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with p* and p* being final state center-of-mass (CMS) mo- 
menta 1 32 1. Neglecting Coulomb and isospin breaking cor- 
rections, which are assumed to be at most a few per cent 1 39 ], 
charge symmetry requires equal strength for the transitions 
n~d — > nn and n + d — > pp. 

Combining optical theorem, charge invariance, detailed 
balance and inserting the parametrisation of the pp — > n + d 
cross section ©, the imaginary part of the n~d —> nn scat- 
tering length is related to the threshold parameter a by 



V>a nD = (! + —)• V>a n - d - 
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The factor (1 + l/S') corrects for the non true absorption 
channels. 

For comparison of pionic-atom and pion-production data, 
the channels pp — > n + d and np — > n°d can be used be- 
cause in the limit of charge independence 2 ■ cr np ^ K a d = 
°~pp^ n+d- Restricting to s waves, in both reactions the same 
transition of the nucleon pair 3 S i(I = 0) — > 3 P\{I = 1) oc- 
curs which corresponds also to pion absorption at rest on 
the deuteron's isospin nucleon-nucleon pair. 



3 Atomic cascade 

After pion capture in hydrogen isotopes, a quantum cas- 
cade starts from main quantum numbers at about n « 16. 
The upper and medium part of the de-excitation cascade is 
dominated by collisional processes (Stark mixing, external 
Auger effect, Coulomb de-excitation). In the lower part X- 
radiation becomes more and more important (Fig.Q} [40 
I4T1I42I. Besides radiative de-excitation by X-ray emission 
all cascade processes depend on the environment, e. g., the 
density of the target gas. 

Stark mixing essentially determines the X-ray yields 
in exotic hydrogen. Because these systems are electrically 
neutral and small on the atomic scale, they penetrate sur- 
rounding atoms and, thus, experience a strong Coulomb 
field. Non-vanishing matrix elements (nlm \ E | nl'm') in 
the presence of the electric field mix atomic states of the 
same principle quantum number n according to the selec- 
tion rules Al = ±1 and Am = 114311 . An induced s state, 
from where pions disapper by nuclear reactions, leads to a 
depletion of the cascade. The Stark mixing rate is propor- 
tional to the number of collisions during the exotic atom's 
life time and explains the strong decrease of the X-ray 
yields with increasing target density BUll . 

In the case of Coulomb de-excitation, the energy re- 
lease of the de-excitation step is converted into kinetic en- 
ergy of the collision partners [44 1. For lower-lying tran- 
sitions a significant energy gain occurs which leads to a 
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Doppler broadening of subsequent X-ray transitions. Doppler 
broadening was directly observed first in the time-of-flight 
spectra of monoenergetic neutrons from the charge exchange 
reaction at rest n~p — > jfin 04511461 . 

A dedicated measurement of the pH(3 p—ls) transition 
has been performed within this series of experiments |47|. 
In pH, where no hadronic broadening occurs, only the Doppler 
broadening exclusively from acceleration effects during the 
atomic cascade contributes and was directly observed as a 
significantly increased line width. The data were used - 
among others - to confirm that the methods used to quan- 
tify the corrections to the line width in pionic hydrogen and 
deuterium are sufficiently well defined [3]. Details may be 
found elsewhere Il48l . 



The acceleration due to Coulomb de-excitation is coun- 
teracted by elastic and inelastic scattering. This leads to a 
rather complex velocity distribution with tails below the 
peaks caused by the Coulomb transitions which are at 12, 
20, 38, and 81 eV for the An = 1 transitions (7 - 6), (6 - 5), 
(5 - 4), and (4 - 3), respectively. 

Cascade calculations have been extended to follow the 
development of the velocity during the de-excitation cas- 
cade and, therefore, predict kinetic energy distributions at 
the time of X-ray emission from a specific level (extended 
standard cascade model ESCM |42|). Calculations exist 
at present only for muonic and pionic hydrogen. Figure|2] 
shows such a prediction for ttH where energies have been 
scaled to the nD case. It turned out in the analysis of the 
pH(3p -Is) line shape that the first ESCM predictions 
could not describe the measured spectrum Il47ll . Therefore, 
a model independent approach was used to extract the rel- 
ative strength of Doppler contributions directly from the 
measured line shape; a method applied successfully first in 
the neutron time-of-flight analysis [46|. Here, the kinetic 
energy distribution is modeled by boxes around the peaks 
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Fig. 1. De-excitation cascade in pionic deuterium. The hadronic 
shift is defined to €\ s = E„ p - Eq ED , i. e., a negative sign as is the 
case in nD corresponds to a repulsive interaction. 
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Fig. 2. Prediction for the kinetic energy distribution of the ttH 
atom at the instant of (3p - Is) transition scaled to the kinemat- 
ics of the 7rD system for a density of 10 bar equivalent. For the 
analysis, the distribution was tentatively approximated by a box 
describing the low-energy component and a second one for the 
high-energy component around 80 eV. The low-energy compo- 
nent reaching 0.55 dominates the distribution (note the broken 
vertical scale). 



generated by Coulomb transitions. One approach as used 
in the analysis of this experiment is included in Figure[2] 

In the case of exotic hydrogen, Auger emission is only 
possible during the collision with other target atoms (ex- 
ternal Auger effect). Auger transitions prefer de-excitation 
steps as small as possible given that the energy gain ex- 
ceeds the binding energy of the electron. Consequently, 
Auger emission contributes mainly in the range n » 6 - 10 
and An = 1 transitions are favored [4042|. Due to the mass 
difference, Doppler broadening from recoil is negligibly 
small compared to the one caused by Coulomb transitions. 

It is known from muon-catalysed fusion experiments 
that during pH + H2 collisions metastable hybrid molecules 
are formed like [(ppp)p]ee [49 50 1 . An analogue process 
is expected in 7rD + D2 collisions, Such complex' are as- 
sumed to stabilise non radiatively by Auger emission. How- 
ever, radiative decay out of molecular states has been dis- 
cussed, and its probability is predicted to increase with nu- 
clear mass 149 50 51 ,52]. Possible X-ray transitions from 
molecular states would falsify the value for the hadronic 
shift determined from the measured X-ray energy. As molec- 
ular formation is collision induced, the fraction of molecules 
formed and the corresponding X-ray rate should depend 
on the target density. Therefore, a possiblity to identify 
such radiative contributions is a measurement of the X- 
ray energy at three different densities (3.3, 10, and 17.5 bar 
equivalent density in this experiment). 



4 Experiment 

The experiment was performed at the 7rE5 channel of the 
proton accelerator at the Paul Scherrer Institut (PSI), Villi- 
gen, Switzerland, which provides a low-energy pion beam 
with intensities of up to a few 10 8 /s (Pig.[3j. Pions of 112 
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Fig. 3. Setup of the nD experiment in the 7rE5 area at the Paul 
Scherrer Insitut (PSI). The roof of the concrete shielding is omit- 
ted tp show the vacuum system connecting the cyclotron trap (up- 
per right), crystal chamber (upper left) and the cryostat of the X- 
ray detector (bottom left). 



MeV/c were injected into the cyclotron trap II HI and de- 
celerated using a set of degraders optimized to the num- 
ber of pion stops in a cylindrical cryogenic target cell of 
22 cm length and 5 cm in diameter in the center of the trap. 
The cell was filled with deuterium gas cooled by means 
of a cold finger. About 0.5% of the incoming pions per 
bar equivalent pressure are stopped in D2 gas. X-radiation 
could exit the target cell axially through a 5 yum thick mylar® 
window. The window foil is supported by horizontal alu- 
minum bars. 

The Johann-type Bragg spectrometer was equipped with 
a spherically bent Si crystal cut parallel to the (111) plane 
having a radius of curvature of R — 2982.2 + 0.3 mm. The 
reflecting area of the crystal having 10 cm in diameter was 
restricted by a circular aperture to 95 mm in diameter to 
avoid edge effects and to 60 mm horizontally to keep the 
Johann broadening small I53,;54l. 

Such a spectrometer is able to measure simultaneously 
an energy interval according to the width of the X-ray 
source when using a correspondingly extended X-ray de- 
tector. Being a pixel device, CCDs are ideally suited for X- 
rays in the few keV range because they combine an intrin- 
sic position resolution with the good energy resolution of 
semiconductor detectors. In addition, the granularity of the 
CCDs allows for efficient background rejection by means 
of pattern recognition (cluster analysis). Photo electrons 
from few-keV X-ray conversion are stopped within a few 
micrometer and, therefore, deposit charge in one or two 
pixels with a common boundary. Beam induced background, 
mainly high energy photons from neutrons produced in 
nuclear pion absorption and captured by surrounding nu- 
clei, leads to larger structures. Together with the massive 
concrete shielding (Fig. [3} such events are efficiently sup- 
pressed (Fig. [5]). 

In this experiment, a 3 x 2 array of charge-coupled de- 
vices (CCDs) was used covering in total 72 mm in height 



and 48 mm in width [55]. Monte Carlo studies show that 
about 2/3 of the intensity of the reflection is covered by 
the height of the array. The detector surface was oriented 
perpendicular to the direction crystal - detector. 

The CCDs' depletion depth of about 30/vm yields an 
optimum quantum efficiency in the 3-4 keV range. A pixel 
size of 40 /mi provides a two-dimensional position resolu- 
tion sufficient to measure precisely the shape of the diffrac- 
tion image. The relative orientation of the individual CCDs 
as well as the exact pixel size was obtained by means of 
an optical measurement using a nanometric grid [56 1. The 
CCDs are operated at a temperature of — 100°C. 

X-ray energies of the ground-state transitions in muonic 
and pionic hydrogen isotopes are in the few keV range, 
where no narrow y lines of sufficient intensity are available 
to be used to determine precisely the function of a low- 
energy crystal spectrometer. However, from highly stripped 
atoms produced in an electron-cyclotron resonance trap 
(ECRIT) narrow X-ray lines are emitted at high rate. 

The spectrometer response was measured using the nar- 
row Ml X-ray line from helium-like argon of 3. 104 keV as 
described in [57 58] yielding a resolution of 436 + 3meV 
(FWHM) when scaled to the energy of the jrD(3p - Is) 
transition of 3.075 keV. This value is close to the intrin- 
sic resolution of 403 + 3 meV given by the rocking curve 
width as calculated for an ideal flat crystal with the code 
XOP 11591. 

The response is constructed from the intrinsic proper- 
ties of the Si (1 1 1) material as calculated from the dynami- 
cal diffraction theory with the code XOP [59|. The rocking 
curve is then convoluted with the experiment geometry by 
means of a Monte-Carlo ray-tracing code (Fig.|4j». Devi- 
ations from this ideal response are introduced by convo- 
luting with an additional Gaussian contribution which was 
determined from the measurement of the Ml transition in 
helium-like argon. 

In the Johann set-up the energy calibration must be 
provided by a reference line of known energy. In order to 
avoid any mechanical changes, only X-rays having about 
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Fig. 4. Transformation of the intrinsic resolution of the Si (111) 
reflection from a flat crystal by folding with the experiment 
geometry and an additional Gaussian describing an imperfect 
mounting. 
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the same Bragg angles are suitable for precision measure- 
ments. The best angular matching for the 7rD(3 p - Is) mea- 
surement was found for the gallium Ka fluorescence lines. 
For energy calibration, the Ga Ka2 transition (9224.484 + 
0.027 eV [60]) was chosen because of the smaller experi- 
mental uncertainty. The energy of the ttD(3/;> -Is) line is 
obtained by the angular distance to the Ga reference line to 
be calculated from the position difference of the reflections 
on the detector and the distance crystal - detector. 

The fluorescence target itself was made of a 25 x20 mm 2 
GaAs plate mounted in the rear part of the gas cell but out- 
side the pion stop volume at 82 mm off the center of the 
cyclotron trap and away from the crystal. The X-rays were 
excited by means of an X-ray tube mounted on a window 
of the cyclotron trap chamber below the gas cell. The value 
for the tabulated Ga Ka energies [60] was obtained from 
measurement also using the compound GaAs lloTl . There- 
fore, a possible chemical shift must not be considered. 

The distance center of the crystal to center of the cy- 
clotron trap was 2100 mm being about 10% outside the 
Rowland circle fulfilling the focal condition R ■ sin 0b- The 
advantage placing the X-ray source off the focal position is 
that one averages over non uniformities of the target. Both 
GaAs and D2 target are large enough that no cuts in the tails 
of the reflection occur. The distance crystal - detector, cho- 
sen close to the assumed ttD focal length, was determined 
to ycD = 1918.1 + 0.5 mm by a survey measurement. 

Alternating measurements of the Ga fluorescence radi- 
ation and the nD line were performed at least once per day 
(Fig-0- The Ga fluorescence X-rays together with two in- 
clinometers mounted at the crystal and detector chambers 
were used to monitor the stability of the line position. De- 
tails on the experimental setup maybe found elsewhere [62 

m. 



5 Analysis 

The X-ray detector raw data are made available as the dig- 
italised charge contents and a position index of the pixel. 
At first, the cluster analysis is performed. As expected, at 
3.1 keV only single («75%) or two pixel events (a 25%) 
contribute. Hot and defect pixels are masked by software. 
The cluster analysed charge (ADC) spectra of the CCDs 
show a pronounced peak originating from 7rD X-rays on a 
largely suppressed background (Pig.[5j. For each CCD, an 
individual energy calibration was performed because of the 
different gain and noise behaviour. The energy resolution 
in terms of charge is determined by means of a Gaussian 
fit. 

Applying an energy cut in the ADC spectra, an addi- 
tional background reduction is achieved. Several ADC cuts 
with a width of 1 cr to Act of the respective Gaussian width 
were used to study the influence of peak-to-background on 
the result for the hadronic broadening. The minimum rel- 
ative statistical error is achieved for a cut of 2.5cr of the 
detector's (charge) resolution. With this energy cut, from 
fits to the position spectra (similar to the one shown in Fig- 
ure |7]- bottom) count rates of 1448+49, 4010+74, and 
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Fig. 5. Background suppression achieved by the cluster analysis 
of the ttD raw data. 

4877+80 were obtained for the nD(3p - Is) line at the 
equivalent densities of 3.3, 10, and 17.5 bar, respectively. 

In the case of the Ga calibration line, compact clusters 
up to size 9 were accepted. Because of the larger penetra- 
tion depth (5 and 105 pm for 3.075 and 9.224 keV in sil- 
icon [64 1, respectively), X-rays also convert at the bound- 
ary of the depletion region where charge diffusion is al- 
ready significant. As the calibration measurements were 
performed without pion beam, the number of background 
events is negligibly small and, hence, no suppression algo- 
rithm is necessary to clean the Ga spectra. 

The resulting hit pattern shows curved reflections on 
the CCD surface. The curvature is corrected for by means 
of a parabola fit (Fig.|6]l before projecting onto the axis of 
dispersion, which is equivalent to an energy axis. 




x/channels 

Fig. 6. 7rD(3p - Is) reflection after correction for curvature. The 
original curvature is indicated by the parabola in the scatter- 
plot of the already corrected data. One channel (pixel) in the di- 
rection of dispersion (x) corresponds in first order reflection to 
76.402+0.001 meV. 



5.1 Energy calibration 

The position of the Ga Ka2 calibration line was determined 
applying a single Voigt profile in the fit. The same proce- 
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dure was used for the tabulated values of the Ga Ka X-ray 
energies [61 ]. 

The ttD line was modeled both with a Voigt profile, 
where Doppler broadening and response function are ap- 
proximated by a single Gaussian, and secondly, the true 
response as determined from the ECRIT data convoluted 
with the imaging properties and including the Doppler con- 
tributions from Coulomb de-excitation (Fig.|7]i. Both meth- 
ods yield the same position value within a few hundreds of 
a pixel. 

A significant correction to the Bragg angle is the in- 
dex of refraction shift |59| because the nD and the Ga Ka 
line is measured in first and third order, respectively. In ad- 
dition, bending and penetration correction due to the sub- 
stantially different energy must be considered. Mechani- 
cal shifts from temperature and vibrations are corrected by 
means of the inclinometer data. 

In total, a systematic error of + 8 meV was deduced 
for the measurement of the nD(3p - 1 s) energy. The com- 
bined statistical error obtained as weighted average from 
the data taken at the three different D2 densities amounts 
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Fig. 7. Ga Ka 2 calibration line and 7rD(3p - Is) transition mea- 
sured at an equivalent density of 10 bar. The narrow structure 
shown inside the nD line represents the resolution of the crystal 
spectrometer. 



to ± 1 1 meV. The largest contribution to the error stems 
from the uncertainty of the energy of the Ga Ko-2 line of 
+ 27 meV. For the transition energy, we find 



= 3075.583 + 0.030 eV. 



(8) 



5.2 Line width 



To the line shape contribute the response of the spectrome- 
ter, the Doppler broadening from Coulomb transitions, and 
the natural line width of the nD(3p -Is) transition. The 
natural line width is given by the width of the 1 s state of 
about 1 eV because nuclear reactions from the 3p level are 
negligible and the radiative width is 28 //eV. 

The model free approach to identify possible Doppler 
contributions uses narrow boxes of e few eV width around 
the kinetic energies as suggested by the energy release of 
Coulomb transitions. Following the experience from the 
analysis of the pH(3p -Is) line shape, one tries to identify 
consecutively individual contributions starting with only 
one contribution of the lowest possible energy, which cor- 
responds to 7tD systems not accelerated or moderated down 
by collisions to energies of a few eV or below. It is suf- 
ficient to approximate the various kinetic energy compo- 
nents by box-like distributions as shown in the analysis of 
the yuH experiment. 

A x 2 analysis, done by means of the MINUIT pack- 
age |65), shows that a low-energy contribution is manda- 
tory. It was found that the kinetic energy must not exceed 
8 eV; a result found independently for taken at 10 bar and 
the 17.5 bar equivalent density (Fig. [8}. The result for the 
natural line width 7Y, turned out to be insensitive to the 
upper boundary of the low-energy box for values < 8 eV. 
Therefore, the low-energy component was fixed to a width 
of 2 eV in the further analysis. 

Noteworthy to mention that the ESCM prediction for 
the kinetic energy distribution for nH{3p - Is) case scaled 
to 7iD kinematics (Fig.[2]i is unable to reproduce the nD(3p- 
1 s) line shape because of rather strong contributions of (4— 
3) and (5-3) Coulomb transitions. Searches for any higher 
energy contributions failed even when using the sum spec- 
trum of the two measurements at 10 and 17.5 bar. This re- 
sult is rather surprising because also the first analyses of 
the new pionic hydrogen data show that at least one high- 
energy component must be assumed to describe the line 
shape of the ttH (4- 1 ), (3 - 1 ), and (2 - 1 ) transitions 11311661 . 
Though a significantly better description of the pH(3 p—l s) 
line shape could be achieved after a recalculation of colli- 
sion cross sections [47 67 68 69 1, there is no explanation 
for the absense of any high-energy components in ttD. 

It has been studied in detail by Monte Carlo which 
magnitude of high-energy energy components maybe missed 
with the statistics achieved in this experiment. The proba- 
bility to miss a high-energy component around 80 e V (box 2) 
corresponding to the (4 - 3) Coulomb transition is dis- 
played in Figure|9] For each set of conditions 400 simu- 
lations were performed. It can be seen, that a high-energy 
contribution of 25% or larger can hardly be missed. For 
10% relative intensity, the chance is about 15% which cor- 
responds to about lcr with respect to the full probability 
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Fig. 8. Search for evidence and width of the low-energy compo- 
nent in the kinetic energy distribution. Diamonds (left x 1 scale) 
are due to the 10 bar data. Squares (right scale) are from a simul- 
tanuous fit to the two spectra taken at 10 and 17.5 bar. 
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Fig. 10. Distributions of relative intensity of a Doppler box 
around Tktn = 80 eV of 10% relative intensity versus the Lorentz 
width r extracted from a set of 400 Monte-Carlo simulations per- 
formed for the statistics of the sum of the 10 and 17.5 bar data. 
The input value f or = 1 1 30 eV. 



distribution allowing a certain variety for the shape of the 
box' extension (Fig.|9). 



2 intensity) and the Lorentz width r reflects the fluctua- 
tions due to the limited statistics (Fig. [Tot. 
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Fig. 9. Top: distributions of relative intensity found from a fit 
with two boxes when a Doppler contribution of 10% and 25% 
relative intensity, respectively, is assumed. The result is based on 
400 Monte-Carlo simulations each performed for the statistics 
of the sum of the 10 and 17.5 bar data. Bottom: probability to 
miss a Doppler contribution at 80 eV as a function of its relative 
intensity. 



Assuming weight zero for box 2 results in an upper 
limit for 7\ s , and using the limit of sensitivity of 10% for 
a second box yields a lower bound for r\ s (- A r slJS ) ac- 
cording to the above-mentioned lcr criterion (Fig. [5]). The 
distribution of the results for the weight a second box (box 



6 Preliminary results 

For extracting the hadronic shift, the pure electromagnetic 
transition energy has been recalculated to be Eqed — 
3077.9062 + 0.0079 eV. The error of this calculation is 
dominated by the uncertainty of the charged pion mass 
(0.0077 eV) GO) where the accuracy of deuteron (0.0012eV) 
and pion radius (0.0010 eV) contribute only marginal. The 
numerical accuracy of this calculation is assumed to be bet- 
ter than lmeVlfTTII. 

Combining measured and calculated QED transition 
energy, we obtain for the hadronic shift 



e u = -(2323 ± 31)meV. 



(9) 



This result is compared to previous measurements in Fig- 
ureQT] Note that the new result for the electromagnetic 
transition energy differs slightly from the values as used 
in the previous experiments 118111911 . 
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Fig. 11. Experimental information on the hadronic shift in pionic 
deuterium. The error of the individual measurements at 3.3, 10, 
and 17.5 bar and, therefore, also of the weighted average is dom- 
inated by the uncertainty of the Ga Ko-2 energy. The results from 
the two previous experiments |18 19 1 are adjusted according to 
the new values for the pure electromagnetic transition energies. 
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Fig. 12. Experimental information on the hadronic broadening 
in pionic deuterium. Data points of previous measurements are 
taken from fUl and fl9l. 



The hadronic broadening is mainly extracted from the 
sum spectrum of the lObar and 17.5 bar measurement us- 
ing only a low-energy component for the kinetic energy 
distribution because no density dependence is identified 
within the experimental accuracy. The combined result 



of about 2 by a remeasurement using a double flat crystal 
spectrometer in case a significant theoretical progress re- 
quires further improvement. The new value for the shift is 
significantly smaller (»4%) than the results reported from 
the two earlier precision measurements. About 40% of the 
discrepancy is due to the new calculation of the 7rD(3 p— 1 s) 
electromagnetic transition energy. 

The theoretical understanding of AW <-» NNn reac- 
tion is continuously increasing. Within the approach of chi- 
ral perturbation theory a study of pion production includ- 
ing the next-to-leading order (NLO) terms yields a NL0 = 
220 pb [ 84 , 87 ] in good agreement with the pionic-deuterium 
result (Fig.rOT>. At NLO the theoretical uncertainty is es- 
timated to about 30%. It is expected to decrease below 
10% within a few years from forthcoming NNLO calcu- 
lations [86]. Hence, the experimental accuracy of 4.2% for 
the hadronic broadening in ttD reaches already the expected 
final uncertainty of the theoretical calculations. 

Noteworthy, that at the 10% level no components from 
high-energy Coulomb transitions could be identified from 
the fit to the line shape. 



A, = 1171 



+ 23 
- 49 



meV 



(10) 
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